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1 Introduction 

In Newtonian gravity the gravitational field of a bounded mass distribution 
is given as a solution of the Poisson equation. As is known, that solution can 
be expanded in a power series of the inverse radial coordinate in the neigh- 
bourhood of infinity, in such a way that some suitable quantities, referred 
to as Multipole Moments (MM), characterize the specific solution. These 
(MM) quantities are related to the structure of the mass distribution by 
means of integrals defined over the source, and hence the solution provides a 
well defined and physically meaningful description of Newtonian gravity. 

When addressing the description of gravity in the vacuum, we look for 
solutions of the Laplace equation, and the above-mentioned multipole ex- 
pansion allows us to identify the arbitrary constants of the general solution 
in terms of the Multipole Moments (MM) of the source. Thus, the general 
well-behavioured solution of the Laplace equation at infinity (or far enough 
from the source) is a series expansion in powers of 1/r the terms of which rep- 
resent the different contributions of each MM that provide both the physical 
characterisctics of the source and a description of the exterior gravitational 
field. 

Let us stress two relevant features of this description of the gravitational 
solutions: first, because of the linearity of the Laplace equation that series 
can be seen as the sum of exact solutions, and hence the partial sums of that 
series become a new solution of the equation. We refer to these partial sums 
as the Multipole Solutions of Newtonian gravity. 

Secondly, we wish to emphasize that the definition of MM gives them a 
well known physical meaning in terms of certain integrals defined over the 
mass distribution. 

These characteristics allow us to consider the multipole expansion series 
as a perturbative approach to the interpretation of the gravitational field by 
following two criteria: the description of the field is performed with respect 
to a point sufficiently far away from the source that higher powers of the 
inverse radial coordinate can be neglected; alternatively, we truncate the 
multipole expansion series at a suitable Multipole Solution because the mass 
distribution of the source provides infinitesimal higher-order contributions. 

Therefore, MM play a fundamental role in the construction of solutions of 
Newtonian gravity since Multipole Solutions can be understood as corrections 
in the desciption of the gravitational field of the first term of the multipole 
expansion (the Monopole solution). This first term of the succession of partial 



sums represents the gravitational field generated by a bounded spherical 
source. The spherical feature characterizes that specific solution univocally 
and, in fact, if the spherical symmetry condition is imposed on the solutions 
of the Laplace equation, then a unique solution of the series is obtained: the 
gravitational monopole. One way to introduce that constraint is by requiring 
the solution to be invariant under the action of the symmetry group 5*0(3). 
In spherical coordinates {r,6,ip}, the infinitesimal generators of that group 
are: 

J2 = -cosip—T + smLpcoiu^- 
od dip 

■ d , d 

Jl = Sm09— - + COS09COt6'7— . 

o9 op 

The invariance condition over a gravitational potential leads to the 

equations — — = 0, — - = 0, and therefore the Laplace equation reduces to 
op 09 

the following equation 

d 



2' 



dr V dr 



0, (2) 



or, equivalent ly, 

2(j)r + rcprr = , (3) 

whose unique solution (different from the trivial constant solution) is simply 

the gravitational monopole = — , (with the imposed condition a > 0, 

r 
from the demand of positive mass). 

This is the point of departure in the line of enquiry pursued here: is there 

any kind of symmetry that describes gravitational multipole solutions? Might 

it be possible to find some differential equations such as ([3]) that restrict 

the solutions of the Laplace equation to those with the prescribed multipole 

characteristic? The spherical symmetry, in the way described above, becomes 

a universal symmetry in the sense that it does not depend on the equation 

we wish to solve namely the Laplace equation. This is not our aim for the 

symmetries we are seeking, but the claim is the existence of some group 

of transformations that will be a symmetry of the Laplace equation, which 

allows us to choose a specific kind of solution from the general equation. 



In what follows we shall see that it is possible to construct a family of 
vector fields that are the infinitesimal generators of a one-parameter group 
of transformations. Furthermore, this group turns out to be a symmetry of a 
system of equations consisting of the Laplace equation and another equation, 
henceforth called supplementary equation. 

Thus, by means of the symmetry, the resolution of the Laplace equation 
when we are looking for multipole solutions simplifies the system of equations 
to be solved, and the problem proves to be easier, since we only need to 
integrate the supplementary equation (which is a very simple ODE in the 
angular variable alone), choosing those given by the correspondent invariance 
condition for the integration functions. 

This work is organized as follows: in section 2, along a first subsection we 
offer a brief review of gravitational multipole solutions in Newtonian gravity, 
and we recall what multipole moments are and their physical meaning. In 
order to avoid too many complicated formulae in the main text, we have used 
an appendix to explain details of these contents. 

In a second subsection we find some differential equations that are ful- 
filled by the partial sums of the series that describe gravitational multipole 
solutions. These equations represent a family of additional conditions to 
the axially symmetric Laplace equation, acting as a restriction to its general 
solution. 

In section 3 we construct a family of vector fields that play the role of 
infinitesimal generators of some one-parameter group of transformations. It 
is proved that these groups are symmetries of a system of equations consisting 
of the Laplace equation and a specific supplementary equation. In a first step, 
we show the so-called Monopole-Dipole symmetry and then we generalize the 
symmetry to the 2^-pole order moment case, where the notion of generalized 
symmetries must be introduced. 

Section 4 is devoted to calculating the groups of transformations for these 
vector fields that generate the symmetries. Finally, in section 5 we show that 
the calculation of group-invariant solutions of the axially symmetric Laplace 
equation leads to the gravitational multipole solutions of Newtonian gravity. 

A Conclusion section contains some comments about the aims achieved 
and possible future generalizations. 

Finally, two appendices are included: Appendix A contains information 
about Classical Multipole Moments, and Appendix B is devoted to complete 
the calculation of the functions involved in the construction of the symme- 
tries. 



2 Multipole solutions in Newtonian gravity 

2.1 The multipole moments and gravitational multi- 
pole solution 

In Newtonian gravity we need to introduce completely symmetric and trace- 
free tensor fields to define the so-called Classical Multipole Moment^j. As 
can be seen in the appendix A, for the case of axial symmetry only one 
component of each tensor is independent and the exterior gravitational field 
generated by an axially symmetric source of mass density p{z, 9) is given by 
the following potential 

°° M 

^ = -GY.:i^Pr.{cose), (4) 

n=0 

where (r, 6) are the radial and polar coordinates of the exterior point with 
respect to any origin on the axes of symmetry; P„ are Legendre polynomials 
and the constants M„ are the multipole moments (I113P of the source. 

Just as the Poisson equation applies for the description of gravity gen- 
erated by a mass distribution, so does the Laplace equation for the vacuum 
case. Thus, if we consider both cases in our configuration space M, the prob- 
lem to solve could be divided into different domains (M = QpUflv) for each 
region as follows: 

Qp, Q^ being the domains of the region containing all the mass distributions 
and vacuum regions respectively. In the Q^ domain, the general solution 
with axial symmetry and good asymptotic behaviour (i.e. decaying to zero 
at infinity or far from the source) 

1 (9 / ^d(j)\ 1 a / . d(f)\ 

J.2 Q^ \ Qrp I J.2 sin 6 86 \ 86 I 



is given by the series: 



oo 
n=0 



^Thc definition and introduction of Classical Multipole Moments is explained with 
details in the Appendix A 



where a„ are arbitrary constants without any physical meaning. Nevertheless, 
since we must demand continuity of the function at the boundary of regions 
Qp and fly, then these constants a„ are (apart from the constant G) identical 
to the multipole moments Mn and acquire the meaning of those quantities. 

2.2 Some properties of multipole solutions 

We shall now show some properties of the series (jl]). The derivative of with 
respect to the variable r proves to be (let us take G = 1): 



•f' = £(" + ^)^^PM = -U\f: nKPM .K^^ (8) 

n=0 n=0 

y = COS 9 , 
and hence we have the following relation: 

rcPr + 4> = Ri (9) 

oo 
R,=Y^nXnPn{y). (10) 



n=l 



The above expression can be read as a linear first-order differential equation 
for the potential 0; it is clear that since Ri, the right part of the equation 
(IHl), does not depend on Aq, the first term of the series (jll) is a solution of the 
corresponding homogeneous differential equation. 

The derivative of with respect to the variable y is: 

oo oo 

0y = - ^ XndyPniy) = "Ai - ^ A„(9j^P„(y) , (11) 

n=0 n=2 

and hence the following relation holds: 

oo oo 

r0^ + + ycpy = ^ nXnPniy) - ^ vKdyPniv) = R2 (12) 

n=2 

00 

R2 = -J2^ndyPn-l{y) , (13) 



n=2 n=2 

00 



n=2 



by means of the property of the Legendre polynomials nPn{y) + P^^i{y) = 
yP^{y). By again considering expression flT2|) as a linear first-order differ- 
ential equation for the potential 0, we realize that since R2 flT^ does not 
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depend on Ai, the partial sum of order 2 (i.e., the sum of the first two terms) 
of the series (jlj) is a solution of the corresponding homogeneous differential 
equation flT2l) . 

Moreover, since the first term of R2 flT3l) is — A2, we can add to the first 
member of the equation flT2|) the second derivative (multiplied by a suitable 
factor) of (p with respect to the variable y (see equation flT^ ) in order to 
obtain a new differential equation for (p with an independent term of higher 
order on A„. By means of a recursive calculation, and considering successive 
derivatives of with respect to the variable y, 

^ n=i ^ n=i+l k=i ^ ' 

where L^^k are the coefficients of the Legendre polynomial of degree n for the 
power k of its variable ?/, we can obtain the following differential equation 
for any value of A^ > 2: 

N 

r0^ + + y0j, + J^/i„(y)9^0 = i?^+i , (15) 

n=2 

where Rn^\ is a series on A^ of ordecl Aat+i, and hniy) are unknown polyno- 
mials of y that we need to introduce to obtain an independent term -Rat+i of 
that order. In fact, starting from (fTOl) and (fT3l) . Rk (with fc > 1) are series 
on A„ which can be written as follows: 

Rk = ^A„Cfc,n(y) , (16) 

n=k 

that is to say, it can be seen that, for all k > 1, Ck,n{y) = with n < k; 
in particular, by using (fT^ and (TT6|) (note that Lj j = (2i — l)!!/i!) one can 
quickly conclude that if one choses Rk recursively in the following way: 



^Thc order of the series must be understood as the smaUest value of the index n for 
a non-null term XnCk,n of the series. Thus, the sum (|16p starts at the value of that 
order. (Note that Ci,o(y) = and 6*2,1 (y) = since Ci,n{y) = nPn{y) and C2,n{y) = 
-dyPn-iiy))- 



then, the A^ term of each Rk+i will vanish for all A; > 1, and therefore from 
equation (HM the functions hn{y) are given by: 

In appendix B we develop the recurrence relation ( TT71) in order to obtain 
the degree and coefficients of the polynomials hn{y), which is needed for 
the construction of the symmetries and the proof of theorems in following 
sections. 

In conclusion, we have that a solution of the corresponding homogeneous 
equation (TT5l) for any A^ is the partial sum of order A^ + 1 of the series 
dlj). Therefore, the homogeneous part of equation (TT5!) can be considered 
as a supplementary condition that must satisfy a solution of the Laplace 
equation in order to be the gravitational multipole solution of order N, or, 
in other words, can be understood as a condition to truncate the series (jl]). 



3 Multipole symmetries 

Let 

Odd 

^^ = ^{x,u)— + T{x,u)— + a{x,u)— (19) 

be a vector field on an open subset M G X x U, where X = M^, with 
coordinates x = {r,y), is the space representing the independent variables, 
and f/ = M, with coordinate u, that represents the dependent variable. 

Let A,^(x,M*-"^) = be a system of u differential equations defined over 
M. As is known [9], if G is a local group of transformations acting on M and 

pr(")v[A^(2;, «("))] =0 , (20) 

whenever A^{x,u^"^) = 0, for every infinitesimal generator v of G, then G 
is a symmetry group of the system, where the vector field pr^^'^v is the n- 
th prolongation of v defined on the corresponding jet space M*^") C X x 
f/('^) , (whose coordinates represent the independent variables, the dependent 
variable, and the derivatives of the dependent variable up to order n) 

pr(")v = v + ya''(a;,M("))-^ , (21) 

^ OUj 



where the summation being over all multi- indices J = (ji, ■■■,jk), with 1 < 

du 
jk ^ '^, ^ ^ k < n, and the notation uj = — — r is used. 

A scalar function "^{x, u) = on M is an invariant of the vector field v 
if the following equation is fulfilled 

(9\1> d"^ d"^ 

d^ du d^ 
and taking into account that -Dr(^) = = t^t; — I" -rr- and Dyi^^) = = 

du or or 
9\1/ du 5\l/ 

— — — — h — — , (where D denotes total derivative) then equation (122|) can be 
ou oy oy 

written in terms of the function u{x) as follows: 

ou ou 
^{x,u)— + t{x,u)- a{x,u)=0. (23) 

3.1 The Monopole-Dipole symmetry 

Let 

d d d _^. 

v = r— + ?/- U-— (24) 

dr dy du 

be a vector field considered as the infinitesimal generator of some group Gmd- 
Looking for invariants \l/(r, y, u) of this vector field (i.e: r'^r+y'^y — u^u = 0) 
we end up with the following condition for the variable u: 

rur + yUy + u = (25) 

which is merely the homogeneous equation (IT^ satisfied by the gravitational 
potential describing the monopole-dipole solution of the Laplace equation. 

We shall now check whether the one-parameter group generated by the 
infinitesimal generator v (1241) is a symmetry group of the Laplace equation. 
In order to do so, we first calculate the second prolongation (J2T1) of v by 
means of formula (see formula (2.39) in P] for details) 

a-^{x, n^^)) = Dj {a{x, u) - rUr - yUy) + ruj^,. + yuj^y , (26) 

where uj^x = -^ — , and we have the following: 

(-21 d d d d d 

pr^ V = V - 2urT: 2mj,- 3^^^- ^Uyy- Sm^j/^ — • (27) 



If this is applied to the Laplace equation (with axial symmetry) 

Au = r'^Urr + 2rUr — 2yUy + (1 — y^)uyy = , (28) 

we obtain: 

pr^^W [Au] = - [Am] - 2uyy . (29) 

Therefore, we can state the following 
Theorem 1: 

The system of equations Ay(a;, m*^"^) = given by 



Ai(x,m(")) = Au = 

A2{X,U^'''^)=Uyy = 



(30) 



where Ai being the Laplace equation fl28l) and A2 the so called supplemen- 
tary equation, admits a symmetry group whose infinitesimal generator is v. 



Proof: 

Since the second prolongation of v acting on the supplementary equation 
is pr'^^W [A2] = —3uyy, and according to fl29|) we have pr'^^W [Ai] = = 
pri'^)y [A2] whenever A^{x, u^"'^) = 0, we conclude with the proof. D 



3.2 2^-pole-order symmetry 

The homogeneous differential equation from expression flTSl) is the equation 
that must be fulfilled by the gravitational multipole solution up to 2^-pole- 
order. In analogy to the previous case, we now introduce a new vector field 
and claim that the condition satisfied by an invariant function of that vector 
field on M reproduces the homogeneous equation (TT5l) . This kind of vector 
field is as follows: 



d d 
or oy 



+ Y.^n{y)d;i 



I (-) 



As can be seen immediately in the above expression, the coefficient function of 

d 
the —- derivative at v depends not only on x and u but also on the derivatives 

ou 
of u. This leads to a significant generalization of the notion of symmetry 
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group, obtained by relaxing the geometric assumption that as long as the 
coefficient function of the vector field depends on x and u such dependence 
will generate a (local) one-parameter group of transformation acting on the 
underlying space M. Noether was the ffist to recognize that it is possible to 
significantly extend the application of symmetry group methods by including 
derivatives of the dependent variables in the infinitesimal generators of the 
transformational. 

Henceforth, we use the term generalized vector field to refer to that kind 
of vector field (13T|) and we shall denote by a[u\ = a [x, u^"'^ any smooth 
differential function depending on x, u and derivatives of u up to order n 
defined for (x,m(")) G M(") C X x f/("). 

Given a generalized vector field such as (1311) . its infinite prolongation is 
the formally infinite sum: 

d d sr-^ J d ,^^, 

p^v = r— + y— + > a -— , (32) 

or oy ^ ouj 

where a^ is given by fl26|) . As is known, a generalized vector field v is 
a generalized infinitesimal symmetry [9] of a system of v differential equa- 
tions Aj^(x,M^"^) = if and only if prv [Aj,] = for every smooth solution 
M = /(x), in direct analogy with the infinitesimal symmetry criterion given 
previously. 

The generalized vector field v (pTl) has the prolongation 

prv = V + a'^ + a^^ + a'^ ^ + a^^-^ + ■ ■ ■ , (33) 

OUr OUy OUrr OUyy 

where the only relevant components for our case are (note that in the follow- 
ing expressions the prime sign is used to denote derivation of functions hn{y) 
with respect to variable y) 

N 
^1" 



-2Ur-J2hn{y)d^Ur (34) 

n=2 

N 

av = -2uy-Y,{h4y)d'^+'u + hMd'^u) (35) 



n=2 



•^A complete discussion of the curious history of generaUzed symmetries can be found 
ini 
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N 



: -'iUrr -^K{y)dyUrr (36) 

n=2 

N 

avv = -3uyy-J2{hn{y)d^^'n + 2hMd^+\ + h:{y)d^u) (37) 



ra=2 



By applying this prolongation to the Laplace equation (!28|) . we obtain: 
prv [Ai] = - Ai [u] - 2uyy + (38) 

N 

- J2 hn{y) (2r9>,, - 2yd';+^u + r^^^^ + (1 - y^)^^"^^^) + 



n=2 

AT 



n=2 

It is straightforward to see that: 
prv [Ai] = - Ai [u] - 2uyy + (39) 

N 

- J2 hn{y) {D-y [Ai [u]] + n(n + l)d;u + 2ni/9;+iw) + 

n=2 

TV 

- E ^^(?/) (-22/^> + (1 - 2/')25r'«) + ^n(2/)(l - y')dyu 

n=2 

and whenever Ai [u] = (i.e., over all smooth solutions u = f{x) of the 
Laplace equation), we have: 

N 

Trr^^ [Ai [u]] = -2uyy - J^ [fi„(y)9> + n„(|/)9;+iw] , (40) 

n=2 

where 

fi„(y) = n{n + l)K{y)-2yh'M + {l-y^)K{y) (41) 

n„(l/) = 2ny/i„(y) + 2(1 - y^)K{y) . (42) 

From the above expressions (l40l) . and considering expressions (I116p . (I117p . 
for the functions hn{y), it can be seen that: 

nkiy) + Uk-iiy)=0 , 3<A;<iV. (43) 

12 



That is to say, it can be proved that the functions hn{y) defined by (11161) 
satisfy (for k >3): 

k{k + l)hkiy)-2yh',iy) + {l-y^)hl + 2{k-l)yhk-i{y) + 2il~y^)h',^,{y) = 

(44) 

And thus, since ^2{y) = ~2, for every vector field v fl3T]) we have calculated 

from (|10])-(|13D that: 

prv [Ai [u]] = -UNd^+\ , N>2 , (45) 

and therefore, since Iliv = 2NyhN{y) + 2(1 — y'^)h'j^{y) ^ 0, we have that 
prv [Ai] =0 if and only if 

an equation that we shall call sitj^jo/ementary equation and denote by A2 \u\. 

With all these results we can state the following: 
Theorem 2: 

The generalized vector field v fl3T]) is a generalized infinitesimal symmetry 
of the system of equations Aj,(x, m*^")) = given by: 



Ai(a;,u(")) = Am = 
Af^(x,M(")) = 9f+^M = 



{m(\. ..(n)\ - n7V+l.. _n (47) 



Proof: 

We shall show that for every smooth solution u = f{x) of the system we 
have pr\ [A,^] =0. 

From equation 046p . it is obvious that pr\ [Ai [u]] = 0. Now, with respect 
to the supplementary equation we have that 



prv 



iV+l 

(TV) 



A^^M =ay-y, (48) 



N+1 



where a^ ' ' ^ is the coefficient of the prolongation of v defined as follows 



N \ 

J2 hn{y)dl'u - rur -yuy\+ rd^+\r + yd^+' 

n=2 J 



^y---y^£)N+ii „. n:^ !, t^.\an^. „„. „.„. \ , ^njv+i„. , ..aN+2 



(49) 
13 



Accordingly, by developing the total derivative we have: 

^-^Vtv + a ^ 

= -{N + 2)d^^\ ~J2[ k J JL^y^'^'^Md^^'u . (50) 

fc=0 ^ ^ n=2 

The first term at the right-hand side of (I3U]) is proportional to the left part 
of equation Ag [u] and therefore, it will be zero for all solution u of the 
system ( H7|) . Moreover, we can see that all derivatives of u with respect to y 
appearing in the sums of (150!) are higher than A^ + 3, because, as we already 
know, the functions hn{y) are polynomials of degree n — 2 with respect to 
y and hence all derivatives d^^^~^hn{y) vanish for A^ + 1 — A;>n — 2, or 
equivalently n + A; < A^ + 3. 

For these reasons, the derivatives of u with respect to y appearing in fl50l) 
can be written as total derivatives of A2 [m], and therefore, we can finally 
write the prolongation acting on A2 [u] as follows: 



prv 



{N) 



(N) 



A^^MJ = -(iV + 2)[ArM«]J+ (51) 

k=0 ^ ^ n=2 

With this argument, it is clear that flSTj) is zero whenever A2 u = 0, and 

so we conclude the proof. D 



4 The group of transformations 

Given a vector field v, the parametrized maximal integral curve T(e, x) pass- 
ing through a; in M is called the flow generated by v. As is known, the flow is 
exactly the same as a local group action on the manifold M or the so-called 
one-parameter group of transformation. 

Therefore, the vector field v is called the infinitesimal generator of the 
action, since by Taylor's theorem in local coordinates we have that: 

T(e,x) = exp{ev){x) = x + e^{x) + 0{e^), (52) 

where ^ = {C,^,- ■ ■ , C,^) are the coefficients of v, and the notation referred to 
as exponentiation of the vector field is used. The orbits of the one-parameter 
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group action are the maximal integral curves of the vector field and there 
is a one-to-one correspondence between local one-parameter groups of trans- 
formations and their infinitesimal generators. 

4.1 Monopole-Dipole group Gmd 

The computation of the fiow generated by the vector field v flT^ . which is the 
infinitesimal generator of the Monopole-Dipole symmetry, is done by solving 
the following system of ordinary differential equations: 

r = ^(eV) |,=o (53) 

y = jyy)V-o 

-u = — (e u) ,=0 , 
de 

which leads to 

exp{e\){x,u) = (eV, e*"?/, e'^w) , (54) 

Thus, the one-parameter group Gmd generated by v of (fT9|l is given by 
the following transformation: 

{x, u) = (eV, e'y, e"'u) ; (55) 

i.e., this symmetry group is a kind of scaling symmetry on M (1 X x U. 

Since the group Gmd was proved to be a symmetry group, if m = /(x) is 
a solution of the system of differential equations fl30l) . so are the functions: 

u = e~'u = e"7(r, y) = e"7(e~V, e~'y) . (56) 

4.2 2^-pole-order groups Gj^j 

With respect to the generalized vector field v ( 131 p . some comments should 
be made about its related group of transformations. First, its one-parameter 
group can no longer act geometrically on the underlying domain M because 
the coefficients of v depend on derivatives of u, which are also being trans- 
formed. Nor can we define a prolonged group action on any finite jet space 
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M*^"\ since the coefficients of pr'^^W will depend on still higher order deriva- 
tives of u than appear in M*^"-' . As is known, the best way to resolve this is 
to define an action of the group on a space of smooth functions as follows [9] : 

[exp(evQ)/] (x) = u{x, e) , (57) 

where u{x, e) is the solution (provided it exists) to the Cauchy problem of 
the system of evolution equations: 

^=Q(x,n(-)) 

u{x,0) = f{x) , (58) 

Q(x,u^'^^) = Q [u] being the so-called characteristic oi the evolutionary vector 
field wq defined as follows ^, [8j: 

VQ = Q N 1^ , (59) 

and associated with any generalized vector field (see formula (5.7) in [H] 
for details). In our case, the characteristic of the evolutionary vector field 
associated with fl3T]) is 

N 

Q[u\ = -u-^ hn{y)dyU - rur - yuy . (60) 

n=2 

As is known [9J, if P [u] is any differential function, and u{x, e) a smooth 
solution to (ISSj) . then the prolongation of the evolutionary vector field deter- 
mines the infinitesimal change in P under the one-parameter group generated 
by vqc|. If we assume convergence of the entire Taylor series in e of the group 
action, we obtain the Lie series: 

P[exp(evQ)/] = ^ ^(prvQ)"P [/] , (61) 

n=0 

where the power n at the prolongation denotes the sucessive application of 
pr\Q n times. In particular, if P \u\ = u then fl6T]) provides the formal series 
solution to the evolutionary system fl58l) : 

u{x, e) = f{x) + eprvQ(/) + \e\pr^,Qf{f) + 0{^) , (62) 



^i.e., from the definition of the group action (|57p we have that 
P[exp{e^rQ)f] = P [/] + epr^Q{P) [J] + 0{^), since ^P [u] = pr^^Q{P). 

16 



d 
where, by means of the formula (see (5.6) in [9j) prvg = y^DjQ [u] — — , 

we have that 

prMf) = Q N 1/ (63) 

(prvQ)2(/) ^ pr^rQ{Q)\f=l-Q-f2hn{y)D^Q-rDrQ-yDyQ]\f ■ 

This formal series fl62|) is no longer of practical use, since we are forced to 
assume that the solution to the Cauchy problem fl58l) is uniquely determined 
provided the initial data f{x) is chosen in some appropriate space of func- 
tions. Accordingly, the resulting flow will be on the above function space. 
Verification of this hypothesis involves overcoming a very difficult problem 
regarding the existence and uniqueness of solutions. And moreover, (162!) 
does not give an explicit expression of the complete transformation group 
but the infinitesimal one, and, as was above mentioned, one has to assume 
its convergence. 

Nevertheless, we have presented here these results because, as we shall 
see in the following section. Lie series solution fl62l) . fl63l) . provides us with an 
argument to introduce a procedure for calculating the G}^ group-invariant 
solutions by considering those solutions of the system of differential equations 
that make the characteristic Q zero. 

An alternative method for calculating the transformation group is now 
approached, as follows [lOj. Let us consider the infinite prolongation fl5^ of 
the generalized vector field v and let us write down the infinite system of 
ordinary differential equations: 

dr dy d uj j 

d e d e d e 

From this system of equations, we can define the fiow of v on the infinite 
jet space to be the solution of the above system with given initial values 
(x, M*^°°^) = (x°, M°, Mj), M° = /(r, y) being a solution of the system of equa- 
tions dlZD: 

ea;p[eprv](x,u(°°)) = (x(e),M(°°)(e)) . (65) 

In our case, from the equation fl64l) we must handle the following system of 
coupled differential equations (we only write the actually relevant equations 
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for our purpose, which is no more than obtaining the flow on the manifold 

M): 



dr 



dy 



de 


de ^ 


du 
de 


N 

-u-y]K{y)u^^ 




n=2 


dUy 

de 


= -{k + l)uf-Df 



N 



E^"(^X 



n=2 



(66) 
(67) 

(68) 



where the notation Uy = dyU has been used. The first two equations in 
(!66|) lead to the already known transformation of the independent variables 

turn out to be 



(' 



r = e'^r^,y = e^y*^), whereas the following ones ([6 



decoupled, as we will see now, when we take into account some considerations 
about the variable u and the functions hn{y)- In order to see this, let us 
consider the above equations 
developed: 



when the total derivative Di has been 



du 



(k 



N 



de 



■(*+iH'-EE 



n=2 j=0 



h^M 



U 



[n+k—j 



(69) 



Let us consider now the derivatives appearing in the right hand side of the 
equation (1^^ . On the one hand, the derivatives of the functions h^iy) (which 
have been denoted by hn (y)) vanish for j > n—2 since the degree with respect 
to y of these functions is n — 2. And on the other hand, the variable u, which 
is a function of the independent variables r and y, is a solution of the system 
of equations fH7|) and therefore, d'^u = 0, Wk > N. Thus, the order of the 



derivatives of u appearing in flB^ fulfills the condition n + k — j < N. With 
these considerations, the non-vanishing terms of the sum of index j in (1691) 
correspond to the following range of the index j: n + k~N<j<n — 2; 
and hence, that equation can be written as follows: 



du 



(k 



N 



n-2 



de 



-{k + r 



u 



{k 



n=2 j=n+k-N 



hliiy) 



u 



(n+k—j 



(70) 



Moreover, for k = N and A; = A^ — 1 in fITOl) the lower value of the 
index j is higher than the upper value and, therefore, the only term in the 



right hand side of the corresponding differential equations of fl70p is the first 
one, whereas for higher values of k, it is straightforward to see that the 
corresponding differential equations can be written as follows: 

duy 



duy 



de 

(N-l 



-(iv + iH 




-Nur^ 






N n-2 


-(iV-a+l)4^-"- 


->;>: 




n=2 j=n-a 




N a-2 , 


-(iV-a+l)Mf-"- 


-\.Y.{ 




n=2 j=0 ^ 



de 



h^r''^\y)v^:'-^ 



N-a 
n — a + ] ^ 

(71) 



where a takes values from 2 to A^ — 1. 

These decoupled equations can be integrated successively with the given 
initial values previously mentioned (x,m*^°°^) = (a;°,M°,Mj), and afterthat we 
can solve, with these solutions, the equation (jHT]) which gives us the group 
transformation of the variable u. 

As a matter of illustration, let us consider for example the case N = 2; 
i.e., let us try to solve the above equations (jSZD, flTT]) in order to obtain the 
action of the group of transformation G\J^ , which represents a symmetry of 
the system of equations fH7|) . The resulting equations for this case are as 
follows: 

du 1 ,^^. 

- = -. + -Uyy (72) 

dUy 

-^ = -2«. (73) 

^ = -3%. , (74) 

and therefore the solution of fl7^ . which is unique derivative involved in fl72|) . 
is given by: 

Uyyix,e)=e-^'uly, (75) 

with the initial condition Uyy{x, 0) = Uyy. And the solution of equation (!72il 
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IS 



uix, e 



''° + l<(i---") 



(76) 



where u{x, e = 0) = u^ has been considered as the initial condition. 

Thus, the one-parameter group generated by v (13T]) in M for the case 
A^ = 2 is given by: 



(x, u) 



e r, e y, e 



«° + 6^.(1 



-2e^ 



(77) 



5 Group-Invariant solutions 

A solution u = f{x) of a system of differential equations is said to be G- 
invariant if it remains unchanged by all group transformations in G, meaning 
that for each element of the group both the function / and that transformed 
by the action of this element agree in their common domains of definition. 
If G is a symmetry group of a system of differential equations A, then we 
can find all the G-invariant solutions to A by solving a reduced system of 
equations, which will involve fewer independent variables than the original 
system. 



5.1 G 



MD' 



■invariant solutions 



According to flS5]) global invariants of the group Gmd are k 
/i = uy, and hence a group-invariant solution takes the form: 



u 



y 



y/r and 



(78) 



Solving for the derivatives of u with respect to r and y in terms of those 
of / with respect to k and so on, and substituing those expressions in the 
Laplace equation, leads to the following linear, constant coefficient differential 
equation: 

/V-2A + 2/ = 0, (79) 

whose general solution is / = an^ + bn. Therefore, the invariant solution of 
the Laplace equation by the action of group Gmd is simply the gravitational 
Monopole-Dipole solution: 



uix] 



a^ + b- 



(80) 
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An alternative way to obtain group-invariant solutions of the system of 
differential equations (H7I) is as follows: First, we solve the supplementary 
equation, which is a PDE but which involves only one coordinate, its general 
solution being u{r, y) = A{r)y + B{r). Now, we make this solution invariant 
under the action of Gmd (l55l) : i.e., the transformed function u(f,y) fl56l) 
functionally becomes the same with u{r,y), and therefore we may conclude 
that the only solutions for A{r) and B{r) are: 

A(r) = 4 , B{r) = -, (81) 

whose transformation by the action of the group proves to be: 

u{r, y) = e-^u{r, y) = e~^ ( ^f^ + ^) = a^ + bl. (82) 

As seen, this method is equivalent to the previous one and proves to be 
easier to develop because we do not need either to search for the charac- 
teristics of the vector field or solve the resulting differential equation (1791) 
after substituing the global invariants and their derivatives in the Laplace 
equation. 

5.2 G]^^ -invariant solutions 

We present now two equivalent methods to obtain the group-invariant solu- 
tions for the case of the generalized vector field v ( 13T1) : firstly, by using the 
explicit expression for the group of transformations obtained by integrating 
fl67|) . flTTl) for any value of A^; and secondly, by using expressions (I62ti63p which 
outline the infinitesimal flow by means of a formal Lie series in terms of the 
characteristic of the evolutionary vector fleld. 

A) For the first method we proceed in analogy with the alternative way 
seen in the previous case of monopole-dipole symmetry. Let us illustrate 
the method by considering, as an example, the case N = 2: we solve the 
supplementary equation of the system of equations ( H71) i.e., u = f{r,y) = 
A{r) + B{r)y + C{r)y'^ and require this solution to be invariant under the 
action of the group GjJ^ ( 1771) . 

Since group G^J was proved to be a symmetry group , if m" = f{x) is a 
solution of the system of differential equations (H7|) . so are the functions: 

1 



u = e "u^ 



f{r,y) + ^fyy{r,y){l-e~'^) 
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(83) 



Thus, the transformed function u{r,y) will be invariant under the action of 
the group if the following equations are fulfilled: 

(84) 

(85) 
(86) 



(87) 



A{r) = e~' 






B{f) = e-^'B{re-') 


C{f) = e-^'Cife-') , 


which have the following unique kind of solutions: 


^W = --^4' B{r) = \. C{r) = ^^ 


The resulting invariant solution turns out to be: 


u{r,y)^ f{r,y) = ' - - ^^ +2/'— +(^)^ 




r r^ 3r-^ 





(88) 

the gravitational multipole solution possessing monopole, dipole and quadrupole 
moments. 

B) The second method to obtain group-invariant solutions has to be with the 
action of the group defined in (157|) . When the formal Lie series ( l62l) . ( l63ll 
was previously outlined, we mentioned that the complete action of the group 
is not nedeed to obtain group-invariant solutions, because an alternative 
procedure arises from those expressions. The argument is supported by the 
following: 

Proposition: 

The functions u = f{x), solutions of the system of differential equations 
A,^(x,M*^")) = given by (H7I) . are group-invariant solutions by the action 
of G]^ (the symmetry group whose infinitesimal generator is given by the 
generalized vector field fIST]) ) if and only if the characteristic of the associated 
evolutionary vector field vq is zero over those functions. 

Proof: 

Since the component of the evolutionary vector field is simply the char- 
acteristic Q[u], if a solution u = f{x) is group- invariant this obviously means 
that vq(u) = 0, and therefore Q{x,u^"'^) = 0. Conversely, from the previous 
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expression fl62ti63p we see that the prolongation of the evolutionary vector and 
all its successive prolongations [prvg)"' depend on Q and its total derivatives, 
and hence the action of the group fl62l) on every solution u = f{x) = u{x, 0) 
that makes Q |„= leads to the equality u{x,e) = u{x,0), i.e. the solution 
is group-invariant because the group transforms the solution into itself. D 



By applying this proposition, we can obtain group-invariant solution of 
the system of differential equations A,^(x, «(")) = (H7I) as follows: We first 
solve Ag and demand that the characteristic of this solution be zero. In 
fact, that condition, Q |„= 0, is exactly the corresponding homogeneous 
equation (TT5l) : 



N 



Q\u=o 



drU - ydyU - J^ K{y)dyU . 



-u — rOrU — 



(89) 



n=2 



The general solution of supplementary equation A2 from flTTj) is: 



(TV) 



N 



u{r,y) = Y^Fk{r)y^ . 



(90) 



/c=0 

If we force that solution to make the characteristic zero, we find that: 

N N 

Y, {F,y^ + Ky'r + kFuy^) + Y. hn{y)d^yU = . 

fc=0 n=2 

It is clear that 



N 



dyU = YP^ 



k\ 



,k—n 



k=n 



{k-n)V 



(91) 



(92) 



and therefore the second sum in the expression (PT]) turns out to be: 



N 



N 



Emi/)5> = E 



n=2 



n=2 



n-2 



i=0 



2^ Fp+n{r) y" 

p=0 ^' 



(93) 



where the following notation has been used for the functions hn{y) 0116p : 



n-2 



^n{y) = Y^^iy' ■ 



(94) 



i=0 
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If we rearrange the sums in powers of the variable y we have: 

N N-2 



Y.K{y)dy = Y.y' 



n=2 



k=0 



' k N 
j=0 p=k+2 



p\ 



{k-j)\ 



(95) 



The equation for the characteristic (19T|) is now: 



N 



N-2 



5^y'=[(A: + l)F, + F^r] + ^ 



y 



k=0 



k=0 



k N 



Y^ J2 F,{r)H, 



-k+j j 



pi 



ik-j)\ 



_j=0 p=k+2 

and so for each power in the variable y we have the following equations: 

A; = A^, A^ - 1 
p\ 



0, 
(96) 



{k + l)Ffc + F^r = 

k N 
j=0 p=k+2 



ik-j)V 



k = 0,---,N-2 



(97) 



whose solutions provide us with the corresponding functions for the potential 
u (!90|) as follows: 



FN{r) 
FN-i{r) 

FJr) 



CN-1 



r 

Ck 



N 



k N 

':p^ ~ ;xw z^ {k-jy. ^ Hp-k+j j J Fi 



{r)r dr 



i=o 
A; = 0, ■ ■ ■ , A^ - 2 . 



p=k+2 



(98) 



(3) 

Let us finish this section with an example. We look for the G\,j group- 
invariant solutions of the Laplace equation. From the above expressions (1981) 
we have that: 

(99) 

(100) 



F.fr) = % 



F,(r) = % 






{H,o + H.,^) j'f^dr 



Cl_ _ 3C3 
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which is the general solution of the equation 2Fi(r) + F[{r)r — -^^ = 0, and 



for /c = in (!98|l . we have: 

Fo{r) = 



6C3 

5H 



^20 / ^ t^r 



--l-s^ (101) 



2 c 
which is the general solution of the equation Fo(r) + Fn(r)r 7 = 0. 

Finally, the invariant solution u{r,y) is: 



Co lC2\ /Ci 3C3\ 2 /C2\ 3 / C3 

£. , £1 , 2 .iM , 2 ,iiM , (102) 

i.e., the partial sum of order 4 of series (jlj) representing the gravitational solu- 
tion constructed with the mass, the dipole, the quadrupole and the octupole 
multipole moments. 



6 Conclusion 

The gravitational multipole solutions of Newtonian gravity are well known, 
and of course it is not the aim of this work to discover them but to introduce 
new insights into the mathematical description of the solutions of the gravi- 
tational potentials with a prescribed number of multipole moments. For the 
case of Newtonian gravity, it is known that these solutions are given by the 
partial sums of the series (j4]). In this work, we have shown that these so- 
lutions can be viewed as group-invariant solutions of certain one-parameter 
groups of transformations Gmd (monopole-dipole case) and G^/ (general 
case). 

The result obtained is a kind of generalization of the fact that the New- 
tonian Monopole is defined by the solution with spherical symmetry of the 
Laplace equation. The existence of some symmetries of the Laplace equation 
is proved, which allows us to restrict the solutions of that equation to those 
with the prescribed MM. In order to do so, it is necessary to append a so- 
called supplementary equation to the Laplace equation in order to provide a 
system that will admit the defined group as a symmetry. 
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The infinitesimal generators of the group as well as the action of the 
group have been constructed for each symmetry (in the general case we cal- 
culate specifically the action of the group G)J to illustrate the procedure). 
And finally it is proved that the group-invariant solutions are exactly the 
gravitational multipole solutions. 

The supplementary equation arises as a constraint imposed on the Laplace 
equation in order to maintain the symmetry of the system of both equations, 
but this restriction itself does not limit the solution of the Laplace equation 
to those desired gravitational multipole solutions if in addition one does not 
demand an asymptotically well-behaved condition at infinity. For example, 
for the case N = 2, the group Gy is a symmetry of the system of equations 
( 147|) whose general solution is: 

^ ^ ai b2y . 2C2P2(y) , , 2 2 p. ^ ..„oN 

u{r, y) =^2- — + ^ + — ^ — + rbiy + -r CiPal^) , (103) 

where the last two terms must be rejected for asymptotic reasons. 

Alternatively, both methods proposed to obtain the group-invariant so- 
lutions for the general case (subsections 5.2A and 5.2B) lead univocally to 
the specific gravitational multipole solutions. In particular, first method pro- 
vides a smart and mathematically standard procedure, with the additional 
advantage that we are dealing with an algebraic condition that is easier to 
solve than a differential onqj. Of course we already know the solution of 
the axially symmetric Laplace equation and the interpretation of the trun- 
cated series from (jl]), but the theoretical result given by the existence of 
some extra symmetries of the Laplace equation is that it provides us with 
a procedure to obtain the gravitational multipole solutions with prescribed 
multipole moments without solving that differential equation. 

A clearly more relevant feature of these results has to be stressed in 
the sense that they serve as a trial for a future generalization to General 



^Note that the general solution above expressed ()103|) requires us to solve, for instance, 
the following coupled differential equations: 

u{r, y) = A{r) + B{r)y + C{r)y'^ 
2rA'{r) + r'^A"{r) + 2C(r) = 
2rB'{r) + r^B"{r) ~ 2B{r) = 
2rC"(r) + r'^C'ir) - 6C(r) = 
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Relativity (GR). Apart from the fact that the non-hnearity of GR gravitation 
does introduce calculation problems, the static and axisymmetric-vacuum 
metrics are described by two metric functions; one of them (which provides 
the other one by means of a quadrature) is actually a solution of the Laplace 
equation. Whereas multipole moments in Newtonian gravity provide physical 
meaning to the coefficients arising from the general solution series ([7]), the GR 
scenario requires complicated computations to obtain the relation between 
those coefficients and Relativistic Multipole Moments (RMM) |1], [1], [2]. 

Some authors have devoted much effort in researching techniques aimed at 
obtaining solutions of stationary axisymmetric vacuum solutions of Einstein 
field equations with prescribed multipole moments. In [5], [1] the Monopole- 
Quadrupole solution was obtained and [5] addresses the MJQ approximate 
solution. In ^ the authors obtained algebraic conditions a^ = an{Mn) re- 
lating the coefficients of the series ([7]) with the RMM, and a metric function 
that is a solution of the Laplace equation is given explicitly. In [6], [7] the 
authors developed a very interesting and useful method for generating the 
coeffcients a„ needed to construct the gravitational multipole solution with 
prescribed RMM. Since the existence and uniqueness of that kind of solu- 
tion in GR can be proved [7], it would be an exceptional goal to establish 
the existence of some kind of "symmetry" that would generalize Birkhoff's 
theorem for a solution with a given set of RMM. The uniqueness theorems of 
partial differential equations are in general very difficult to solve and a topic 
of research in Mathematical Physics, and they would not strictly be the aim 
of future works. However, the proposal outlined here consists of consider- 
ing the algebraic relations a„ = an{Mn) from a different mathematical point 
of view, hopefully searching for the existence of some kind of "symmetry" 
rather than a boundary condition problem from which those relations arose. 

Since it is Laplace's equation that we wish to solve, the problem could be 
oriented towards finding the appropriate system of coordinates or some basis 
of functions where the already known gravitational multipole solutions, as 
in the case of Newtonian gravity, could be considered as the only solutions 
satisfying suitable symmetry conditions. 
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7 Appendix A 

Let us consider the newtonian gravitational potential of a mass distribution 
with density p{z), given by the following solution of the Poisson equation 

Hx) = -gJ ^p{z)d'z , (104) 

where G is the gravitational constant; the integral is extended to the volume 
of the source; z is the vector that gives the position of a generic point inside 
the source, and R is the distance between that point and any exterior point 
P defined by its position vector x. Let us now make an expansion of this 
potential in a power series of the inverse of the distance from the origin to 

the point P (r = |s|) by means of a Taylor expansioiij of the term — around 

the origin of coordinates, where R = \/{x'^ — z^){xi — Zi). The result is the 
following: 

$(f) = -^-Gf;i^g--n.,...n, ; n^ ^ ^ , (105) 

1=1 

where M represents the total mass of the source, i.e., 

M= p{z)£z , (106) 

Jv 

and the quantities Q*i'"*' are completely symmetric and trace- free tensor 
fields defined as follows: 



Qn-M ^ ^21 - l)\\T f z'' 
Jv 



..z''p{z)d^z , (107) 



where the symbol T denotes the subtraction of the traces. 

By using spherical harmonic functions it is possible to give a better rep- 
resentation of the multipole moments from the above-mentioned expansion 

of the function — if it is written in terms of the Legendre polynomials as 

R 
follows: 

1 °° 1 

R = E;JTT^'^Kcos/5) , (108) 

1=0 



^see [TTj for details 
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where z ^\z\ and /3 is the angle defined between the position vector of the 
exterior point x and the position vector of the interior point of integration z. 
Let us write the Legendre polynomials in terms of the spherical harmon- 
ics: 

. m=-l 

A-7T 

Piicosp) = ^^^ Yl Yrie,<i>)Yr{eA) , (109) 

m=—l 

where {6, 0) and {6, 0) are the angular spherical coordinates of the position 
vectors x and i* respectively. Thus, the gravitational potential $(x) generated 
by a mass distribution outside the source is given by 

rU oo m=l 

^i^) = -^-^-Gj2^f^,;jr.J:DrYr{0A) , (110) 

1=1 m=—l 

D™ being some quantities defined as follows: 

Dr= [ z%*"'{ej)p{z)d^z . (Ill) 

Jv 

For each value of /, there exist 2/ + 1 complex quantities Dp. Nevertheless, 
since V;~™ = (— 1)™Y;*™, only 2/ + 1 real quantities are independent. 

From expressions flllOj) and (11051) it is clear that the completely symmetric 
and trace- free tensors Q*^'"*' are equivalent to the 2/ + 1 real independent 
quantities l| D 

^Q^-''n,, . . .n, = ^ Yl Dryr((^A) • (112) 

m=—l 

Finally, if we consider axial symmetry the number of quantities needed to 
define the multipole moments decreases. The only spherical harmonic that 
provides a non-vanishing integral in the expressions of D™ is the correspond- 
ing yj°, since the solution does not depend on the azimuthal angle. Therefore, 
only one quantity remains to define the multipole moment of order /: 

Mi = 2n [ [ z^+^p{z,e)Pi{cose)smededz , (113) 



m=l 



^A completely symmetric tensor Ti^,,,i^ of rank r, in a manifold of dimension 3 only 
has (r + 2)(r + l)/2 independent components. In addition, the trace-free condition, i.e., 
Ti-^...ir9^^^^ = imposes r{r — l)/2 constraints, and therefore only 2r + 1 components of a 
completely symmetric and trace- free tensor in dimension 3 are independent. 
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p{z^ 9) being the density of the axially symmetric distribution; the integral is 
extended to the volume of the source, and therefore z represents the radius 
of the integration point and 9 the corresponding polar angle. 

8 Appendix B 

The functions hn{y) of equation (TT^ have been defined by (TT51) in terms of 
some quantities Ck,n{y) from flT6l) . By using formula flT^ in the recurrence 
relation flT7|) we immediately get 

C^.+i,n(2/) = C,,„(y) - ^^;i^9jP„(y) , n>k . (114) 

This recurrence relation is solved by just plugging it into itself as many times 
as needed, but nevertheless it is not necessary to solve all the quantities of 
these double-index functions Ck,^ since we are only looking for the quantities 
Cn^n- As has been said, by developing the recurrence relation (I114p (with 
increasing value of k) starting from C2,n we easily see that 

CkAy) = C.M - E (2^'-f)!! ^^^"(^^ , n>fc . (115) 

Finally, by taking k = n and considering (fT8l) and C2,„ = —dyPn-i{y) we get 
the following relation to calculate the functions hn{y)'- 

n-l 

dyPn-i + Y,hk{y)d'yPn{y) 



Kiy) 



i2n - 1) 



fc=2 



^ > 3, h2{y) = --. 



h2{y) 



(116) 
From these relations (11161) we can see that the functions hn{y) are poly- 
nomials in the variable y of degree n — 2: 

1 

~3!! 

2 

h,{y) = -^(i + V) 

hM = |y(32/ + %^) 

h,{y) = -^{l + l2y^ + %y') (117) 

30 



Acknowledgments 

This work was partially supported by the Spanish Ministry of Education and 
Science under Research Project No. FIS 2006-05319. This work would not 
have been possible without the ideas provided by Professor Jesus Martin, 
whom I wish to thank for helpful discussions spread out over a good number 
of years. I also wish to thank Dr. Javier Villaroel for his comments about the 
bibliography and Dr. Alberto Alonso and Dr. Miguel Angel Gonzalez-Leon 
for very fruitful insights. 

References 

[1] Geroch, R. (1970) J. Math. Phys., 11, 2580. 

[2] Hansen, R.O. (1974) J. Math. Phys. 15, 46. 

[3] Hernandez- Past or a, J.L., Martin, J. (1993) Class. Quantum Grav. 10, 
2581. 

[4] Hernandez- Past or a, J.L., Martin, J. (1994) Gen. Rel. and Grav.. 26, 

877. 

[5] Hernandez- Past or a. (2006) Gen. Rel. and Grav. 38, 871. 

[6] Backdahl, T., Herberthson, M., (2005) Glass. Quantum Grav. 22, 3585 
Herberthson, M., (2004) Glass. Quantum Grav. 21, 5121 
Backdahl, T., Herberthson, M., (2006) Glass. Quantum Grav. 23, 5997 

[7] Backdahl, T., Herberthson, M., (2005) Glass. Quantum Grav. 22, 1607 

[8] Gaeta, G., (2005) Int. J. Geom. Meths. Phys. 2 1081 Glass. Quantum 
Grav. 17, 3617. 

[9] Olver, P. J., (1993) Applications of Lie Groups to Differential equations. 
Graduate Texts in Mathematics. Springer- Verlag Ed. New York. 

[10] Anderson, R.L., Ibragimov, N.H. (1979) Lie-Bdcklund transformations 
in Applications SIAM Studies in Appl. Math., 1, Philadelphia. 

[11] Hernandez-Pastora, J.L., (1996) Ph.D. Relativistic gravitational fields 
close to Schwarzschild solution. Universidad de Salamanca. 



31 



